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Comment on "Why is the DNA denaturation 
transition first order?" 

Recently, Kafri, Mukamel and Peliti (KMP) 1] ex- 
tended the classical Poland-Scheraga (PS) model for the 
denaturation transition of DNA jj. According to PS, 
the energetic binding of bases in the double helix com- 
petes with the entropic contribution of denatured loops, 
implying that the nature of the transition depends on the 
exponent c in the statistical weight il.{2k) ^ s^k~^ for a 
closed loop of length 2k: for 1 < c < 2 it is of second 
order, while for c > 2 it is of first order Fisher, 
taking the effects of self-avoidance within a denatured 
loop into account, found c = dv ^ 1.766 in d = 3, i.e., 
a second order transition 0]. In Ref. it was obtained 
that the exponent c is modified if additional effects of 
self- avoidance between a denatured loop and the vicinal 
double helices are included. For a single loop within two 
strands of double helix, it was found that Q 


c — dv — 2(T3 


2.115 


d = 3 


(1) 


i.e., a first order transition, where (T3 is a topological 
exponent related to a 3-vertex of a polymer network 

This conclusion is valid in the asymptotic scaling limit 
for long flexible, self-avoiding chains, i.e., each of the 
three segments going out from a vertex must be much 
longer than the persistence length £p of this segment 
(even though individual values of the ip may be differ- 
ent). If this condition is fulfilled, the analysis for the 
PS-inspired model in Ref. 1] is consistent How- 
ever, we point out in this Comment that it does not apply 
to the chains typically used in experiments 7J : the DNA 
double helix being quite rigid, we expect the transition 
in such systems to be of second order. 

The typical length of DNA used in experiments varies 
from about 100 to 5000 base pairs (bp), the latter cor- 
responding to a whole viral DNA Q. In such DNA, a 
chain "monomer" m, which corresponds to a bead in 
a freely jointed chain, represents one persistence length 
£p. For the single strand in a denatured loop, typically 
£p(L) ~ 40 A (roughly 8 bases), whereas for the double 
helix £p(H) - 500A (100 bp) Even if one assumes 


that a segment of 10 monomers to(H) of the double he- 
lix is already long enough to be sufficiently close to the 
asymptotic scaling limit for long chains (which is hope- 
lessly optimistic |9( , and also much less than taken in the 
simulations |^), one can at best place 5 loops even on 
the longest chains such that the flexibility condition is 
not violated. However, with a maximum number of only 
5 loops, the system is governed by finite size effects, and 
the analysis in Ref. 1] is no longer valid. 

Conversely, only 80 bases (40 bp) are needed to form 
10 monomers m(L) such that a denatured loop can be 
considered as sufficiently flexible. The longer chains in 
the experiments can thus exhibit a fairly large number of 
such loops, if the segments of the double-stranded helix 
between them are allowed to be of the order of £p(H) (see 
below). It is therefore justifled to neglect the entropy of 
the double-stranded helix 0|. Following this picture, a 
vertex with three outgoing legs would not tie together 
three flexible chains, but rather two flexible chains (be- 
longing to the loop) and one rigid rod (the double helix). 
However, a rigid rod represents an irrelevant object for 
flexible, self-avoiding chains in the scaling limit for d = 3 
[Toj , which implies for the present case that (T3 in Eq. ^ 
is replaced by zero. This, in fact, reproduces the original 
value c ~ 1.766 Q, i.e., a second order transition. 

Finally, we note that at the denaturation transition of 
real DNA the average length ^ of bound helical segments 
between loops is of the order of a few hundred bp, thus 
comparable with ^p(H) and much larger than ^p(L) (see, 
e.g., Eq. (9.129) in Ref. 0). This large value of ^ reflects 
the fact that the denaturation transition of real DNA 
is highly cooperative; in the PS model, this is generally 
modeled by the (non-universal) cooperativity parameter 
fo ^ 1 13) H uM which enters the statistical weight of 
loops in the partition function, leading to ^ ~ I/cto 3> 1 
(see Eqs. (9.1), (9.129) and (9.115a) in Ref. 
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